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Abstract 


There  is  currently  a  desire  to  use  Chebyshev  polynomials  to  fit  terrain 
elevation  data.  Such  a  fit  would  create  a  surface  function  that  exactly  fits  the 
known  elevations,  and  would  describe  an  elevation  at  any  point  on  that 
surface.  This  note  questions  the  appropriateness  of  using  Chebyshev 
polynomials  for  this  purpose,  as  opposed  to  linear  interpolation  or  use  of  a 
cubic  spline.  A  set  of  elevations  in  one  direction  is  used  to  illustrate  a  point 
that  large  transitions  in  elevation  influence  the  coefficients  in  the  polynomial 
fit  and  contribute  spectral  energy  to  points  far  from  the  transition  area.  It 
argues  that  a  linear  interpolation  process,  or  a  cubic  spline  interpolation, 
would  be  more  appropriate. 
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1.  Introduction 


There  is  a  concern  whether  a  Chebyshev  polynomial  curve  fit  to  elevation 
data  is  an  appropriate  method  for  interpolating  the  elevation  values  for 
points  between  those  that  are  known,  such  as  those  given  by  the  Defense 
Mapping  Agency  Digital  Terrain  Elevation  Data,  level  1  coverage.  In  brief, 
a  Chebyshev  pol)momial  is  defined  as  (Hamming,  1986;  Press  et  al,  1994) 

T„(x)  =  cos(«  cos~\x)) 
or 

W  =  1 

Tjfx)  =  a: 

T2ix)  =  2x^-1 

T^ix)  =  4:^-3x 

T„(x)  =  2xT„_j  (x) - T„_2  (x)  n>2  . 


A  curve  fit  using  Chebyshev  polynomials  takes  on  the  form 


where  the  Chebyshev  polynomial  coefficients,  c*,  are  determined  from 

C](  =  -x7^  f{x^)  1J._2  (Xf^)  , 

W=1 


where 

N  =  total  number  of  data  points  to  fit, 

f(x)  =  exact  function  described  by  a  Chebyshev  polynomial  curve  fit, 
f(x„)  =  data  point  at  location  n,  and 

c*  =  Chebyshev  polynomial  coefficient  for  polynomial  k-1. 

By  taking  advantage  of  the  polynomial's  orthogonality,  we  can  determine 
the  coefficients  in  a  simpler  manner  by  substituting  for  x„  with 


which  results  in 


cos 


N  J 
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This  last  relationship  is  important  because  it  gives  the  Chebyshev  pol)mo- 
mial  a  Fourier  transform  quality.  In  essence,  the  coefficients  are  the  product 
of  a  discrete  Fourier  transform  of  the  elevation  data.  This  Avill  later  give 
rise  to  problems  of  anti-aliasing  when  the  curve  fit  is  created. 


2.  Use  of  a  Chebyshev  Polynomial  Curve  Fit 

The  rationale  for  using  a  Chebyshev  polynojnial  to  fit  a  curve  is  that  it  will 
exactly  fit  all  points  considered  in  the  fit.  This  type  of  fit  uses  the  spectral 
information  of  the  known  points  to  determine  the  behavior  of  a  function  in 
between  and  over  all  known  points.  The  Chebyshev  polynomial  is  also 
very  convenient  to  use  since  the  polynomial  fit  at  an  arbitrary  point  need 
not  be  determined  from  the  entire  set  of  coefficients.  There  is  a  point  in  the 
summation  where  additional  terms  are  negligible.  Often  times,  only  about 
20  percent  of  the  terms  are  needed  to  converge.  In  the  terrain  data  set,  it  is 
quite  easy  to  produce  a  matrix  of  coefficients  that  would  represent  a 
Chebyshev  "surface."  An  elevation  value  would  be  simply  calculated  from 
this  array  of  coefficients.  Unfortunately,  there  are  some  problems. 

Terrain  data  need  not  be  continuous.  Cliffs,  overhangs,  river  gorges,  can¬ 
yons,  mesas,  or  essentially  any  rapid  transition  in  elevation  along  an 
isolattitude  or  isolongitude  would  present  a  problem  in  the  Chebyshev  fit. 
Remember  that  the  Chebyshev  coefficients  are  very  similar  to  the  coeffi¬ 
cients  generated  by  a  Fourier  transform.  Consider  a  set  of  discrete  points 
with  somewhat  abrupt  changes  in  value,  such  as  depicted  in  figure  1. 
Straight  lines  are  drawn  between  the  points  to  illustrate  the  effect  of  a  lin¬ 
ear  interpolation.  As  is  easily  seen,  there  are  some  large  transitions,  with 
the  largest  occurring  between  points  1  and  2.  A  rough  estimate  of  the  slope 
is  somewhat  steep  at  Ay/Aoc  =  -0.46.  These  points  are  not  an  unreasonable 
expectation  of  hilly  terrain  leading  to  a  valley.  In  fact,  they  would  represent 
a  likely  cross-section  of  an  avenue  for  ground  troop  movement  or 
encampment. 


Figure  1.  Elevation 
data  along  an 
arbitrary 
isolattitude. 


2 


Horizontal  distance  (m) 


A  Chebyshev  polynomial  fit  results  in  a  curve  that  is  shown  in  figure  2. 
One  can  immediately  see  the  effect  of  anti-aliasing  in  the  first  1000  m  of  the 
terrain  profile,  which  is  brought  on  by  the  highly  transient  regions  contrib¬ 
uting  to  the  value  (or  energy)  of  the  high  wave  number  terms.  These  terms 
are  later  reflected  in  the  oscillations  near  the  base  of  the  hill  as  continuing, 
yet  decaying,  in  amplitude.  In  fact,  everywhere  in  this  section  of  terrain, 
the  Chebyshev  fit  has  introduced  high  spectral  content  in  the  fit.  This  is  the 
occurrence  in  the  Fourier  transform  of  a  discontinuous  function.  This  is 
typically  not  true  in  the  termination  of  hills  into  a  valley. 


3.  Cubic  Spline  Interpolation 

So  the  question  becomes,  is  it  still  more  reasonable  to  fit  terrain  elevation 
with  a  Chebyshev  polynomial  rather  than  use  a  linear  interpolation,  or  is 
there  something  better?  Perhaps  the  use  of  a  cubic  spline  evaluated  from  a 
few  surrounding  points  would  be  a  better  way  to  accoimt  for  local  trends, 
still  gain  a  nonlinear  estimate  of  the  terrain  elevation,  and  yet  not  have  ef¬ 
fects  from  terrain  that  is  quite  distant  from  the  area  of  interest.  The  cubic 
spline  presents  a  compromise  between  the  linear  interpolation  and  a 
Chebyshev  pol)momial  fit.  The  changes  between  points  are  represented  by 
a  rovmded  transition  to  better  represent  hills  and  bluffs.  Although  precipi¬ 
tous  drops  are  still  smoothed  and  remain  underrepresented  in  the  terrain 
fit,  the  results  are  probably  better  than  a  simple  linear  interpolation.  Look¬ 
ing  at  the  same  terrain  data,  but  using  a  cubic  spline  to  fit,  we  see  in  figure 
3  that  these  two  criteria  are  met;  namely,  only  local  effects  are  considered, 
and  the  curvature  (in  this  case,  based  on  the  calculated  second  derivative) 
is  used  to  estimate  the  deviation  from  linearity. 

A  cubic  spline  interpolation  is  probably  better  for  estimating  elevations, 
especially  in  the  case  of  multiple  interpolations  between  actual  data  points 
(desired  resolution  greater  than  actual  grid  spacing).  A  cubic  spline  would 
resolve  a  smooth  transition  about  the  data  point.  If  elevation  across  several 
points  is  needed  (desired  resolution  less  than  actual  grid  spacing),  linear 
interpolation  is  most  likely  adequate,  since  fine  terrain  features  are  filtered 
out  by  the  larger  grid  spacing,  and  average  values  representing  the  pre- 


Figure  2.  Chebyshev 
polynomial  fit  to 
terrain  data. 
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Figure  3.  Cubic  spline 
fit  to  terrain  data. 
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dominant  terrain  of  the  area  are  more  reasonable  for  a  fluid  dynamic  com¬ 
putation  boundary.  The  computational  cost  of  evaluating  the  Chebyshev 
polynomial  surface,  and  then  extracting  the  interpolated  values,  is  much 
more  costly  than  performing  a  cubic  spline  interpolation,  as  the  cubic 
spline  interpolation  is  more  costly  than  linear  interpolation. 


4.  Conclusion 

By  this  assessment,  the  cubic  spline  interpolation  is  beneficial  only  if  grid 
spacing  less  than  twice  the  elevation  data  grid  spacing  is  needed.  Error 
associated  with  linear  interpolation  is  probably  of  the  same  order  as  the  er¬ 
ror  associated  with  a  cubic  spline  for  larger  grid  spacing.  A  Chebyshev 
polynomial  fit  near  a  high  relief  area  would  result  in  excessive  elevation 
oscillations  that  may  not  be  realistic.  It  is  recommended  that  a  cubic  spline 
interpolation  process  be  used  for  fine  grid  interpolation  and  that  linear  in¬ 
terpolation  be  used  for  coarse  grid  interpolation. 


References 

1.  R.  W.  Hamming,  Numerical  Methods  for  Scientist  and  Engineers,  2*’'^  edition, 
Dover  Publications  Inc.,  New  York  (1986). 

2.  W.  H.  Press,  S.  A.  Teukolsky,  W.  T.  Vetterling,  and  B.  P.  Flannery,  Numerical 
Recipes  in  FORTRAN.  The  Art  of  Scientific  Computing,  2”^*  edition,  Cam¬ 
bridge  University  Press  (1994). 


4 


Distribution 


Admnstr 

Defns  Techl  Info  Ctr 
Attn  DTIC-OCP 

8725  John  J  Kingman  Rd  Ste  0944 
FTBelvoirVA  22060-6218 

Central  Intellgnc  Agency 

Dir  DB  Standard 

Attn  GE  47  QB  (10  copies) 

Washington  DC  20505 

Defns  Intllgnc  Acgy 
Attn  DT  2  Wpns  &  Sys  Div 
Washington  DC  20301 

Dir  of  Defns  Rsrch  &  Engrg 
Attn  DD  TWP 
Washintgon  DC  20301 

Commanding  Officer 
Attn  NMCB23 
6205  Stuart  Rd  Ste  101 
FT  Belvoir  VA  22060-5275 

DARPA 
Attn  Techl  Lib 
3701  N  Fairfax  Dr 
Arlington  VA  22203-1714 

Dir  of  Chem  &  Nuc  Ops  DA  DCSOPS 
Attn  Techl  Lib  (2  copies) 

Washington  DC  20310 

Dpty  Assist  Scy  for  Rsrch  &  Techl 
AttnSARD-TR  R  ChaitRm3E476 
AttnSARD-TT  K  Kominos 
AttnSARD-TT  F  Milton  Rm  3E479 
AttnSARD-TT  C  NashRm3E479 
The  Pentagon 
Washington  DC  20310-0103 

Hdqtrs  Dept  of  the  Army 
AttnDAMO-FDQ  D  Schmidt 
AttnDAMO-FDQ  MAJM  McGonagle 
400  Army  Pentagon 
Washington  DC  20310-0460 


US  Army  Engrg  Div 
AttnHNDEDFD 
PO  Box  1500 
Huntsville  AL  35807 

US  Army  Mis  &  Spc  Intllgnc  Ctr 

Attn  AIAMS  YDL 

Redstone  Arsenal  AL  35898-5500 

US  Army  NGIC 
Attn  Rsrch  &  Data  Branch 
220  7thStretNE 
Charlottesville  VA  22901-5396 

US  Army  Nuc  &  Cheml  Agency 
7150  Heller  Loop  Ste  101  (3  copies) 
Springfield  VA  22150-3198 

US  Army  Rsrch  Lab 
Attn  SLCRO-D 
PO  Box  12211 

Research  Triangle  Park  NC  27709-221 1 

US  Army  Strtgc  Defns  Cmnd 
Attn  CSSD  H  MPL  Techl  Lib 
AttnCSSDHXM  Dr  Davies 
PO  Box  1500 
Huntsville  AL  35807 

Chief  of  Nav  OPS  Dept  of  the  Navy 
Attn  OP  03EG 
Attn  OP  985F 
Washington  DC  20350 

Nav  Civil  Engrg  Lab 
Attn  Code  L51  J  Tancreto 
Port  Hueneme  CA  93043-5003 

Officer  in  Charge 
Nav  Construction  Battlion  Ctr 
Civil  Engrg  Lab 
Attn  Code  L31  Techl  Lib 
Port  Hueneme  CA  93041 

Nav  Rsrch  Lab 

Attn  Code  2027  Techl  Lib 

Washington  DC  20375 


Distribution 


Ofc  of  Nav  Rsrch 

Attn  Code  23  A  Faulstick  (2  copies) 
800  N  Quincy  Stret 
Arlington  VA  22217 

US  Nav  Academy 
Attn  Tech  Lib 
AttnP  McCoy 
Attn  LTD  Fuller 
572  Holloway  Rd 
Annapolis  MD  21402-5002 

US  Nav  Academy 
AttnR  Smith 

363  Rickover  Hall  590  Holloway  Rd 
Annapolis  MD  21402-5002 

MIT 

Attn  Tech  Lib 
Cambridge  MA  02139 

Univ  of  Maryland 
Attn  Engrg  &  Techl  Lib 
AttnJ  Duncan 
AttnJ  Wallace 
Attn  K  Kreiger 


Univ  of  Maryland  (cont’d) 

AttnL  Ong 
AttnP  Bernard 

Rm  2168  Engrg  Classroom  Bldg 
College  Park  MD  20742-5121 

US  Army  Rsrch  Lab 

Attn  AMSRL-CI-LL  Tech  Lib  (3  copies) 
Attn  AMSRL-CS-AL-TA  MaU  &  Records 
Mgmt 

Attn  AMSRL-CS-AL-TP  Techl  Pub 
(3  copies) 

Attn  AMSRL-IS-EE  A  Wetmore 
Attn  AMSRL-IS-EE  D  Garvey 
Attn  AMSRL-IS-EE  J  Martin 
Attn  AMSRL-IS-EE  K  Deacon 
Attn  AMSRL-IS-EE  R  Cianco 
Attn  AMSRL-IS-EE  R  Loucks  (5  copies) 
Attn  AMSRL-IS-EE  R  Meyers 
Attn  AMSRL-IS-EE  R  Pinnick 
Attn  AMSRL-IS-EE  S  Niles 
Attn  AMSRL-IS-ES  COL  Price 
Attn  AMSRL-IS-ES  P  Gillespie 
Adelphi  MD  20783-1197 


6 


REPORT  DOCUMENTATION  PAGE 


Form  Approved 
0MB  No.  0704-0168 


Public  reporting  burden  for  this  collection  of  information  is  estimated  to  average  1  hour  per  response,  induding  the  time  for  reviewing  instructions,  searching  existing  data  sources, 
gathering  and  maintaining  the  data  needed,  and  completing  and  revievring  the  collection  of  information.  Send  cornments  regarding  this  burden  estimate  or  any  other  aspect  of  this 
collection  of  information,  including  suggestions  for  reducing  this  burden,  to  Washington  Headquarters  Services,  Directorate  for  Information  Operations  and  Reporte,  121 5  Jefferson 
Davis  Highway.  Suite  1204,  Arlington.  VA  22202-4302.  and  to  the  Office  of  Management  and  Budget.  Papenwork  Reduction  Project  (07&4-0188),  Washington,  DC  20503. 


1.  AGENCY  USE  ONLY  (UavB  blank) 


i  REPORT  DATE 

December  1 996 


3.  REPOBTTVPE  AND  DATES  COVERED 

Interim,  from  1  to  31  Oct.  1996 


4.TTrLE  AND  SUBTITLE 

5.  FUNDING  NUMBERS 

Chebyshev  Polynomial  Fit  for  Terrain  Elevation 

DA  PR:  B53A 

PE:P61120 

6.  AUTHOR(S) 

Richard  B.  Loucks 

7.  RERFORMIN6  ORGANIZATION  NAME(S)  AND  ADDRESS<ES) 

U.S.  Army  Research  Laboratory 

Attn:  AMSRL-IS-EE 

2800  Powder  Mill  Road 

Adelphi,  MD  20783-1197 

8.  PERFORMING  ORGANIZATION 

REPORT  NUMBER 

ARL-TN-86 

A.  SPONSORINGMONITORING  AGENCY  NAME(S)  AND  ADDRESS(ES) 

U.S.  Army  Research  Laboratory 

2800  Powder  Mill  Road 

Adelphi,  MD  20783-1197 

10.  SPONSORING/MONITORING 

AGENCY  REPORT  NUMBER 

11.  SUPPLEMENTARY  NOTES 

AMS  code:611102.53A11 

ARL  PR:  7FEJ60 

12l.  DISTRIBUTION/AVAILABILITY  STATEMENT 

12b.  DISTRIBUTION  CODE 

Approved  for  public  release;  distribution  unlimited. 

13.  ABSTRACT  (Maximum  ^  vforxis} 

There  is  currently  a  desire  to  use  Chebyshev  polynomials  to  fit  terrain  elevation  data.  Such  a  fit  would 
create  a  surface  Unction  that  exactly  fits  the  known  elevations,  and  would  describe  an  elevation  at 
any  point  on  that  surface.  This  note  questions  the  appropriateness  of  using  Chebyshev  polynomials  for 
this  purpose,  as  opposed  to  linear  interpolation  or  use  of  a  cubic  spline.  A  set  of  elevations  in  one 
direction  is  usecf  to  illustrate  a  point  that  large  transitions  in  elevation  influence  the  coefficients  in  the 
polynomial  fit  and  contribute  spectral  energy  to  points  far  from  the  transition  area.  It  argues  that  a 
linear  interpolation  process,  or  a  cubic  spline  interpolation,  would  be  more  appropriate. 


14.  SUBJECT  TERMS 


Chebyshev,  terrain,  fit 


IS.  NUMBER  OF  PAGES 

13 


17.  SECURfTY  CLASSIFICATION 
OF  REPORT 

Unclassified 


18.  SECURITY  CLASSIFICATION 
OF  THIS  PAGE 


Unclassified 


19.  SECURITY  CLASSIFICATION 
OF  ABSTRACT 

Unclassified 


20.  UMITAT10N  OF  ABSTRACT 

UL 


NSN  7540^1-280-5500 


Standard  Form  298  (Rev.  2-89) 
Prescribed  by  ANSI  Std.  Z3d<l8 
298-102 


7 


